Compact Central WENO Schemes for Multidimensional Conservation Laws

We present a new third-order central scheme for approximating solutions of systems of conservation laws in one and two space dimensions. In the spirit of Godunov-type schemes,our method is based on reconstructing a piecewise-polynomial interpolant from cell-averages which is then advanced exactly in time. In the reconstruction step, we introduce a new third-order as a convex combination of interpolants based on different stencils. The heart of the matter is that one of these interpolants is taken as an arbitrary quadratic polynomial and the weights of the convex combination are set as to obtain third-order accuracy in smooth regions. The embedded mechanism in the WENO-like schemes guarantees that in regions with discontinuities or large gradients, there is an automatic switch to a one-sided second-order reconstruction, which prevents the creation of spurious oscillations. In the one-dimensional case, our new third order scheme is based on an extremely compact point stencil. Analogous compactness is retained in more space dimensions. The accuracy, robustness and high-resolution properties of our scheme are demonstrated in a variety of one and two dimensional problems.Comment: 24 pages, 5 figure

Quinpi: Integrating Conservation Laws with CWENO Implicit Methods

Many interesting applications of hyperbolic systems of equations are stiff, and require the time step to satisfy restrictive stability conditions. One way to avoid small time steps is to use implicit time integration. Implicit integration is quite straightforward for first-order schemes. High order schemes instead also need to control spurious oscillations, which requires limiting in space and time also in the linear case. We propose a framework to simplify considerably the application of high order non-oscillatory schemes through the introduction of a low order implicit predictor, which is used both to set up the nonlinear weights of a standard high order space reconstruction, and to achieve limiting in time. In this preliminary work, we concentrate on the case of a third-order scheme, based on diagonally implicit Runge Kutta (DIRK) integration in time and central weighted essentially non-oscillatory (CWENO) reconstruction in space. The numerical tests involve linear and nonlinear scalar conservation laws

Central Schemes for Nonconservative Hyperbolic Systems

In this work we present a new approach to the construction of high order finite volume central schemes on staggered grids for general hyperbolic systems, including those not admitting a conservation form. The method is based on finite volume space discretization on staggered cells, central Runge-Kutta time discretization, and integration over a family of paths, associated to the system itself, for the generalization of the method to nonconservative systems. Applications to the one and the two layers shallow water models as prototypes of systems of balance laws and systems with source terms and nonconservative products respectively, will be illustrated

One- and Multi-dimensional CWENOZ Reconstructions for Implementing Boundary Conditions Without Ghost Cells

We address the issue of point value reconstructions from cell averages in the context of third-order finite volume schemes, focusing in particular on the cells close to the boundaries of the domain. In fact, most techniques in the literature rely on the creation of ghost cells outside the boundary and on some form of extrapolation from the inside that, taking into account the boundary conditions, fills the ghost cells with appropriate values, so that a standard reconstruction can be applied also in the boundary cells. In Naumann et al. (Appl. Math. Comput. 325: 252-270. https://doi.org110.1016/j.amc.2017.12.041, 2018), motivated by the difficulty of choosing appropriate boundary conditions at the internal nodes of a network, a different technique was explored that avoids the use of ghost cells, but instead employs for the boundary cells a different stencil, biased towards the interior of the domain. In this paper, extending that approach, which does not make use of ghost cells, we propose a more accurate reconstruction for the one-dimensional case and a two-dimensional one for Cartesian grids. In several numerical tests, we compare the novel reconstruction with the standard approach using ghost cells

ARBITRARY ORDER FINITE VOLUME WELL-BALANCED SCHEMES FOR THE EULER EQUATIONS WITH GRAVITY

This work presents arbitrary high order well balanced finite volume schemes for the Euler equations with a prescribed gravitational field. It is assumed that the desired equilibrium solution is known, and we construct a scheme which is exactly well balanced for that particular equilibrium. The scheme is based on high order reconstructions of the fluctuations from equilibrium of density, velocity, and pressure, and on a well-balanced integration of the source terms, while no assumptions are needed on the numerical flux, beside consistency. This technique also allows one to construct well-balanced methods for a class of moving equilibria. Several numerical tests demonstrate the performance of the scheme on different scenarios, from equilibrium solutions to nonsteady problems involving shocks. The numerical tests are carried out with methods up to fifth order in one dimension, and third order accuracy in two dimensions

An all speed second order well-balanced IMEX relaxation scheme for the Euler equations with gravity

We present an implicit-explicit well-balanced finite volume scheme for the Euler equations with a gravitational source term which is able to deal also with low Mach flows. To visualize the different scales we use the non-dimensionalized equations on which we apply a pressure splitting and a Suliciu relaxation. On the resulting model, we apply a splitting of the flux into a linear implicit and an non-linear explicit part that leads to a scale independent time-step. The explicit step consists of a Godunov type method based on an approximative Riemann solver where the source term is included in the flux formulation. We develop the method for a first order scheme and give an extension to second order. Both schemes are designed to be well-balanced, preserve the positivity of density and internal energy and have a scale independent diffusion. We give the low Mach limit equations for well-prepared data and show that the scheme is asymptotic preserving. These properties are numerically validated by various test cases

Efficient asymptotic preserving schemes for BGK and ES-BGK models on cartesian grids

This work is devoted to the study of complex flows where hydrodynamic and rarefied regimes coexist. This kind of flows are found in vacuum pumps or hypersonic re-entries of space vehicles where the distance between gas molecules is so large that their microscopic behaviour differ from the average behaviour of the flow and has be taken into account. We then consider two models of the Boltzmann equation viable for such flows: the BGK model dans the ES-BGK model. We first devise a new wall boundary condition ensuring a smooth transition of the solution from the rarefied regime to the hydrodynamic regime. We then describe how this boundary condition (and boundary conditions in general) can be enforced with second order accuracy on an immersed body on Cartesian grids preserving the asymptotic limit towards compressible Euler equations. We exploit the ability of Cartesian grids to massive parallel computations (HPC) to drastically reduce the computational time which is an issue for kinetic models. A new approach considering local velocity grids is then presented showing important gain on the computational time (up to 80$\%$). 3D simulations are also presented showing the efficiency of the methods. Finally, solid particle transport in a rarefied flow is studied. The kinetic model is coupled with a Vlasov-type equation modeling the passive particle transport solved with a method based on remeshing processes. As application, we investigate the realistic test case of the pollution of optical devices carried by satellites due to incompletely burned particles coming from the altitude control thrusters